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ABSTRACT 


The performance of two novel interferometer antenna architectures for high- 
resolution, wideband direction finding are investigated. The first configuration incorporates 
a Symmetrical Number System (SNS) encoding of the interferometer amplitude response 
(symmetrical folding waveform). The second configuration incorporates a Residue Number 
System (RNS) encoding of the interferometer phase response (saw-tooth waveform). The 
residue architectures serve as a source for resolution enhancement in an interferometer array 
by decomposing the analog spatial filtering operation into a number of parallel sub¬ 
operations (moduli) that are of smaller computational complexity. Each sub-operation only 
requires a precision in accordance with the size of the modulus. A much higher resolution 
is achieved after the N moduli are used and the results of these low precision sub¬ 
operations are recombined. A four-element, 3 channel array using the moduli set mj = 3, 
m 2 = 4 and m 3 = 5 was constructed in a ground plane using rectangular waveguide 
elements with a center fi'equency of 8.5 GHz. Experimental results are compared with the 
simulation results to demonstrate the advantages of this approach. The frequency response 
of the RNS array is investigated numerically. To correct the quantization errors due to any 
frequency offset, a fast correction algorithm is derived and is shown to have excellent 
results over a wide bandwidth. 
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1. INTRODUCTION 


A. OVERVIEW OF DIRECTION FINDING (DF) SYSTEMS 

Direction Finding (DF) systems provide an emitter’s bearing, or angle of arrival that 
can be used as an invariant sorting parameter in the deinterleaving of agile pulsed radar 
signals and in separating closely spaced communications emitters. In addition, the 
conservation of jamming power in power-managed Electronic Attack (EA) systems 
depends on the ability of the associated Electronic Warfare Support (ES) system to measure 
the direction to the victim emitter. Direction-of-arrival measurement systems fall into three 
main categories: scanning beam, simultaneous-multiple-beam, and interferometer 
techniques. The mechanically scanning beam requires only a single antenna. However, this 
technique has the disadvantage of a low probability of intercept. The simultaneous- 
multiple-beam system uses an antenna or several antennas, forming a number of 
simultaneous beams, thereby retaining the high sensitivity of the scanning beam approach 
while providing fast response. However, it requires many parallel receiving channels, each 
with full frequency coverage. Interferometer systems, based on phase-comparison 
techniques, have the advantage of a fast response but generally use wide coverage antennas 
which result in low sensitivity. In addition, they require relatively complex microwave 
circuitry that must maintain a precise phase match over a wide frequency band under 
extreme environmental conditions. When high accuracy is required (on the order of 0.1 to 1 
degree), wide baseline interferometers are used with ambiguity resolving circuitry. [Ref. 1] 

The two primary interferometer techniques used for DF are the amplitude- 
comparison method and the phase-comparison method. Generally the phase-comparison 
method has the advantage of greater accuracy, but the amplitude-comparison method is 
widely used due to its lower complexity and cost. The basic interferometer geometry is 
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illustrated in Figure 1.1 whereby a plane wave arriving at an angle 6 is received by one 
antenna element earlier than the other due to the difference in path length. The time 
difference can be expressed as a phase difference 

A(() = ^-^sin6 (1-1) 

where d is the spacing between antennas, and X is the wavelength. For a two-element 
interferometer, the spacing d must be X/2 or less to provide unambiguous, or single lobe 
±90° coverage. 



Figure 1.1: Phase Interferometer Principle. 
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The spacing between antennas must be established for the highest frequency to be 
received. Interferometers employing multiple antenna elements are called multiple baseline 
interferometers. Usually the receiver consists of a reference antenna and a series of 
companion antennas. The spacing between the reference element and the first companion 
antenna is X,/2; other secondary elements are placed to form pairs separated by 1, 2, 4, and 
8 wavelengths. This multiple baseline interferometer provides a binary angle-of-arrival 
measurement where each bit of the measurement supplies a more accurate estimate of the 
emitter’s angle of arrival. [Ref. 1] 

Another kind of multiple baseline interferometer is the nonharmonic. In 
nonharmonic interferometers, no pair of elements provides completely unambiguous angle- 
of-amval information over the complete field of view. For example, the initial spacing in 
the nonharmonic interferometer might be X, while the next companion element spacing is 3 
X/2. Ambiguities are resolved by tmth tables, and hence the accuracy is set by the spacing 
of the widest baseline antenna pair. Nonharmonic interferometers have been implemented 
over nine-to-one bandwidths (2 to 18 GHz) with rms accuracy from 0.1 to 1 degree and 
without ambiguities over ±90°. The principal advantage of the nonharmonic interferometer 
over the harmonic interferometer is the increased bandwidth for unambiguous coverage 
[Ref. 2]. 

B. PRINCIPAL CONTRIBUTIONS 

This thesis investigates the performance of two new interferometer antenna 
architectures for high-resolution wideband direction finding. Hie first configuration 
incorporates a symmetrical number system (SNS) encoding of the interferometer amplitude 
response (symmetrical folding waveform). Tlie SNS can serve as a source for resolution 
enhancement in an interferometer by decomposing the analog spatial filtering operation into 
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a number of parallel sub-operations (moduli) that are of smaller computational complexity. 
Each sub-operation symmetrically folds the antenna’s amplitude response with folding 
period equal to twice the modulus. Thus, each sub-operation only requires a precision in 
accordance with that modulus. A much higher resolution is achieved after the N different 
SNS moduli are used and the results of these low precision sub-operations are recombined. 
Simulation results are used to numerically evaluate the performance of a 4-element (3 
channel) array using moduli mj=3, m2=4 and m3=5. The antenna was constmcted in a 
ground plane using rectangular waveguide elements with a center frequency of 8.5 GHz. 
Experimental (anechoic chamber) results were obtained and compared with the simulation 
results to demonstrate the feasibility of the approach. The second antenna configuration 
incorporates a residue number system (RNS) encoding of the interferometer phase 
response (sawtooth folding waveform). The phase response folds with folding period 
equal to the modulus. Simulation results indicate that the RNS architecture antenna 
provides a more robust solution to high-resolution DF. The frequency response of the RNS 
array is investigated numerically. Quantization errors in the resolved direction of arrival 
when the intercepted frequency is other than the center frequency are demonstrated. A fast 
correction algorithm that uses the known intercepted frequency (e.g., from an intercept 
receiver) is derived. Tfre performance of the algorithm is quantified and shown to have 
excellent results over a large bandwidth. 

C. THESIS OUTLINE 

Chapter n details the mathematical development of the linear antenna array. A two 
element interferometer is shown as a special case. The array factor associated with the 
linear array is derived. Chapter III, is devoted to the encoding of amplitude response with 
the SNS. First the SNS encoding procedure is briefly introduced. The SNS antenna 
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architecture is then derived with simulation results shown to demonstrate the array 
performance. Construction of the antenna hardware is described and experimental results 
are shown. Chapter IV examines encoding the phase response with the RNS. The RNS 
encoding procedure is briefly described and simulation results are shown. In Chapter V the 
antenna frequency response, antenna size and respective accuracy are investigated for the 
RNS antenna design. The correction algorithm is derived using the intercepted frequency. 
Chapter VI presents a summary, some concluding remarks and recommendations. The 
Appendix provides the Matlab™ code that simulates the output of both the SNS and RNS 
anterma for direction finding, as well as the calculations of the quantization errors. 
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11. LINEAR ARRAY ANTENNA THEORY 


A. OVERVIEW 

Several radiation elements can be arranged in space and interconnected to produce a 
directional radiation pattern. Such a configuration of multiple radiating elements is referred 
to as an array antenna, or simply an array. In an array of identical elements, the radiation 
pattern is determined by the type of individual elements used, their orientations, their 
positions in space, and the amplitude and phase of the current feeding them. To provide 
specialized radiation patterns, it is necessary that the fields from the elements of the array 
interfere constructively (add) in the desired signal directions and interfere destructively 
(cancel each other) in the remaining space. Ideally this can be accomplished, but practically 
it is only approached. Arrays can be found in many geometrical configurations. The most 
elementary is that of a linear array in which the array element centers lie along a straight 
line. The elements may be equally or unequally spaced. When the array element centers are 
located in a plane, it is said to be a planar array. Examples of planar arrays are circular and 
rectangular arrays in which the element centers are disposed on a circle, or are contained 
within a rectangular area, respectively. Another class of arrays is that of conformal arrays. 
In this application the array element locations must conform to some nonplanar surface 
such as that found on an aircraft or missile. [Ref. 3] 

The simplest and one of the most practical arrays is formed by uniformly 
distributing the elements along a line. This arrangement is the conventional linear array 
configuration. Considering that all of the array elements are identical, the principle of 
pattern multiplication applies. This is stated as: “The electrical field pattern of an array 
consisting of similar elements is given by the product of the pattern of one of the elements 
(the element factor (EF)) and the pattern of an array of isotropic point sources with the 
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same locations, relative amplitudes and phases as the original array (the array factor (AF))” 
[Ref. 4]. That is, 


E(e,(i)) = EF(e,(t))-AF(e) (2-1) 

where E(0, is the total radiated (or received) electric field intensity in the direction (6, 

<t)). 

The array factor may be thought of as a sum of isotropic point sources located at the 
center of each array element. The array factor is that factor of the radiation pattern which is 
found from the element currents (amplitudes and phases) and their locations. On the other 
hand, the element pattern is that factor of the radiation pattern determined by the individual 
properties of an element (its current distribution and orientation in space). 

B. LINEAR ARRAY ANTENNA BEAMFORMING 

A typical linear array composed of similar element antennas is shown in Figure 
2.1. The output of each array element may be controlled in amplitude and phase as 
indicated by the phase shifters and attenuators. In addition to the amplitude and phase 
control imposed on each element represented by I„ (n = 0, 1, 2,..., N^-1), there is relative 
phase shift between the waves arriving at the antenna elements due to their positions in 
space and the angle of arrival of the wave. Furthermore, the pattern of each element leads to 
a response that varies with the angle of arrival of the incoming plane wave. If it is the same 
for each element, the principle of pattern multiplication allows us to consider it as one 
multiplicative factor m the total pattern. The remaining angular dependence of the pattern is 
called the array factor, and it is determined only by the element positions and flieir 
amplitudes and phases represented by L. The array in Figure 2.1 is shown as a receiving 
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array; however, the receiving pattern is the same as the transmitting pattern by reciprocity 
[Ref. 4]. 



Figure 2.1: A Typical Linear Array. The Output Currents are Summed in a Summing 

Device Before Entering the Receiver. 

Since the pattern of an array obtained by neglecting die patterns of the individual 
elements of the array is the array factor, then if we replace each element of the array with an 
isotropic point source, the resulting pattern is the array factor. The isotropic point source is 
an idealized antenna occupying a point in space and radiating uniformly in all directions. 
The radiation fields of a point source at the origin of a spherical coordinate system are 
proportional to I^e ^/47cr where is the current of the point source, k is the wave-number 
equal to and r the radius of the sphere [Ref. 4]. The far-field pattern is obtained from 
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the angular dependence of the radiation fields. Thus the pattern of a point source is constant 
and given by 


AF = V (2-2) 

As elements are added to the array, each one with a different current, it is necessary to 
account for their relative field strengths as determined by their element currents. 

The array factor for the array in Figure 2.1 is derived using the array in Figure 2.2 
which has point sources for the array elements in place of the actual elements. 



Figure 2.2: Equivalent Configuration of the Array in Figure 2.1 for Determining the Array 
Factor. The Elements of the Array are Replaced by Isotropic Point Sources. 
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The array factor for this receiving array is then the sum of the point source responses 
{e^'''%e^'''',...,e''‘''‘‘''} each weighted by the amplitude and phase shift {Ig, 
introduced in the transmission line connected to each element. Here are the 

phases of an incoming plane wave at the elements locations designated 0,1,..., -1. 

Usually these phases are relative to the coordinate origin; that is, the phase of the wave 
arriving at nth element leads the phase of the wave arriving at the origin by \|rj^ _j. [Ref. 4]. 
Thus, the array factor of the array shown in Figure 2.2 is given by 

AF = Ioe''^»+I,ej'^'+ (2-3) 

Now, consider a linear array of elements as shown in Figure 2.3. 



Figure 2.3: Equally Spaced Linear Array of Isotropic Point Sources. 


The elements are equally spaced with a distance between elements equal to d. The angle 6 is 
that of an incoming plane wave relative to the axis of the array. The isotropic sources 
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respond equally in all directions, but when their outputs are added together (each weighted 
according to I„) a directional response is obtained. The phase of the wave arriving at die 
origin (n=0) is arbitrarily set to zero, \)ro=0. The incoming wave at element 1 arrives before 
those at the origin since the distance is shorter by an amount dsinS. The corresponding 
phase leads at element 1 relative to those at element 0 is kdsin0. This process continues and 
Equation (2-3) becomes 

N.-l 

AF = lo -f--" = (2-4) 

n=0 

Now consider the array to be transmitting. If the current has a linear phase 
progression (i.e., the relative phase between adjacent elements is the same), we separate 
out this phase explicitly as 


I 


n 



(2-5) 


where the n+7 element leads the nth element in phase by (x. Then Equation (2-4) becomes 


Define \|/ as 



jn(kd sinS+a) 


n=0 


\|/ = kd sin6 + a 


(2-6) 


(2-7) 


then 


N.-l 

AF=2A.e'"». 

n=0 


( 2 - 8 ) 
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A very important special case of equally spaced linear arrays is that of the uniformly excited 
array. This is an array whose element current amplitudes are identical, so that 

Aq — Aj — Aj = ••• = Aj^^_j. (2-9) 

The array factor from Equation (2-8) is then 

N.-l 

AF = Ao =Ao[ l + e^'^+e^^'^+•••-»-(2-10) 


A few short steps are required to sum this geometric series. First multiply the last equation 
by e^'*' to obtain 


AF e^'^ = Ao [e^'^ + e'"'^ + • • • + e^‘ '^]. 


( 2 - 11 ) 


Subtracting Equation (2-11) from Equation (2-10) gives 


or 


AF(l-ej'^) = Ao(l-e^‘'^) 

/ 


AF = Ao 


l_ejN.v 

l-e^'^ 


( 2 - 12 ) 

(2-13) 


This may be rewritten in a more convenient form as follows 


pjN.V/Zr jN,v/2 -jN,V/2'| 

AF = A -_- = A -_ ^ 

0 e'*''-! ® e^'*''^ e^v/2_g-jv/2 > 


so that 
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(2-14) 


AF = A 


_j(N.-l)v/2 

0® 


sin(N^\|//2) 

sin(\|r/2) 


The phase factor of the array is given by If the array was centered about the 

origin, this phase factor would not be present since it represents the phase shift of the array 
phase center relative to the origin. Neglecting the phase factor, Equation (2-14) becomes 


AF = Ao 


sin(N^V}r/2) 

sin(\i//2) 


(2-15) 


The maximum value of the array factor is equal to A^N^ and occurs when \|/=0. To 
normalize it so that the maximum value is equal to unity (if Ao=l) and is independent of the 
number of elements. Equation (2-15) is divided by 


1 sin(N,\|//2) 
sin(\|//2) 


(2-16) 


Now let 6(, be the corresponding value of 0 for which the array factor is maximum. Then 
from Equation (2-7) we have \|/ = kd sinG^ + a = 0, or 


a = -kdsin0o = 


2jt d 

~X~ 


sin 00 . 


(2-17) 


Equation (2-17) represents the element-to-element phase shift in the excitation currents 
necessary to produce an array factor main beam maximum in a direction dg relative to the 
line along which the array elements are aligned. This main beam scanning by phase control 
mode can be implemented by substituting Equation (2-17) into Equation (2-7) obtaining 
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V = kd(sin e - sin 00 ) = (sin 0 - sin 0o). 

A 


(2-18) 


Setting the phases of the array in this way is called “steering the beam” [Ref. 3]. 


C. TWO-ELEMENT INTERFEROMETER ANTENNA 

For a two-element interferometer antenna (N^=2), Equation (2-16) becomes 


AF ■ 

2sin(y/2)' 


Assuming the two-element configuration is centered about the origin of the reference 
system, and using the trigonometric identity 


sin(2v) = 2sin(\)/)cos(\|/), 


Equation (2-19) can be simplified further becoming 


AF„. 


sin(\)//2)cos(\|//2) 

sin(\j//2) 


or 


AF„„^ = cos(\|//2). 

A plot of the array factor using Equation (2-21) is shown in Figure 2.4 for d=6L 


(2-19) 


( 2 - 20 ) 


( 2 - 21 ) 
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Normalized AF Amplitude 


1 



Figure 2.4: Array Factor Amplitude for a Two-Element Linear Array Separated by a 

Distance d=6X,. 


The phase response of the two-element interferometer antenna with d=6X. is derived from 
Equation (2-18) with 0 q=O and -90° < 6 < 90° as illustrated in Figure 2.5a. The sawtooth 
waveform results by constraining \|r between -Jt and tc. 
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Phase Response (rad) 



Figure 2.5a: Phase Response of a Two-Element Linear Array Separated by a 

Distance d=6 X,. 


Figure 2.5b shows the resultant phase response folding waveform that will be encoded 
using the RNS. 
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Figure 2.5b: Phase Response of a Two-Element Linear Array Separated by a Distance 
d=6 X (\|t constrained to vary between -% and 7t). 


Note the sawtooth nature of this folding waveform. Also note the change in folding period 
as the input direction of arrival approaches ±90°. 
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III. SYMMETRICAL NUMBER SYSTEM ENCODING OF THE 

AMPLITUDE RESPONSE 


A. THE SYMMETRICAL NUMBER SYSTEM (SNS) 

The optimum SNS is composed of a number of Pairwise Relatively Prime (PRP) 
moduli mj. The integers within each SNS modulus are derived by amplitude analyzing a 
symmetrically folded waveform with the period of the waveform equal to twice the PRP 
modulus, i.e., 2tr\ [Ref. 5]. Given m, the integer values within twice the individual 
modulus are given by the row vector 


x = [0,l,-,m-l,m-l,-,l,0]. (3-1) 

From (3-1) the required number of threshold levels is mj-l. By analyzing Equation (3-1) it 
is clear that the integers within each modulus do not form a complete system of length 2m, 
by themselves due to the presence of ambiguities. However, the ambiguities that arise 
within the modulus can be resolved by using an N moduli set (mj, m 2 ,...,mj^) and 
combining the values from aU the channels together. By doing this, the SNS becomes a 
complete system having a one-to-one correspondence with the residue number system. For 
N equal to the number of PRP moduli, the dynamic range of this scheme is [Ref. 5,6] 

M = nm,. (3-2) 

i=l 

This dynamic range is also the position of the first repetitive moduli vector and the total 
number of quantization levels without ambiguities. As an example consider the moduli set 
mi=3, m2=4, and m3=5. The first repetitive moduli vector occurs at an input equal to 60, 
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and the number of quantization levels without ambiguities is also 60 (0 to 59). Table 3.1 
illustrates the SNS for the moduli set (3, 4, 5). Figure 3.1 shows the symmetrical folding 
waveforms and SNS states for mj=3, m2=4, and m3=5. The individual states within the 
SNS represent the number of threshold levels crossed by the folding waveform. A total of 
nij-l threshold levels are required for each folding waveform. 


Normalized 

Input 

mi=3, 

Moduli 

m2=4, 

m3=5 

0 

0 

0 

0 

1 

1 

1 

1 

2 

2 

2 

2 

3 

2 

3 

3 • 

4 

1 

3 

4 

5 

0 

2 

4 

6 

0 

1 

3 

7 

1 

0 

2 

8 

2 

0 

1 

9 

2 

1 

0 

10 

1 

2 

0 

11 

0 

3 

1 

12 

0 

3 

2 

« 

52 

1 

3 

2 

53 

0 

2 

3 

54 

0 

1 

4 

55 

1 

0 

4 

56 

2 

0 

3 

57 

2 

1 

2 

58 

1 

2 

1 

59 

0 

3 

0 

60 

b 

3 

0 


Table 3.1: An Example to Illustrate the Optimum SNS 
Using Moduli Set mi=3, mi=4, m3=5. 
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Array Factor Output 


□ ; SNS states (number of threshold levels crossed) 



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

Normalized Input Direction of Arrival 


Figure 3.1: Optimum SNS Waveforms and Output Codes for mj=3, m2=4, and m3=5. 


































B. SNS ANTENNA ARCHITECTURE 


Figure 3.2 shows a schematic representation of the SNS antenna architecture for 
instantaneous direction finding. The array consists of four elements. Element 1 is combined 
with elements 2, 3, and 4. The three pairs of elements (1, 4), (1, 3), and (1, 2) form a 
three-baseline modular interferometer with mi=3, ni2=4 and m3=5 respectively. 



Figure 3.2: Schematic Representation of the Proposed SNS Antenna 
Architecture for Combat Direction Finding. 
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The array factor amplitude response from each pair of elements (corresponding to modulus 
m^) is quantized using mj-l comparators. The comparator threshold levels are adjusted to 
midlevel quantize the folded amplitude response into the SNS format. The SNS-to-decimal 
logic block then recombines each channel and converts the SNS representation to a more 
familiar digital oulput 

The distance between each pair of elements must be determined in order to obtain 
the correct number of folds needed to cover the entire field of view (from 0 to ±90°) for the 
chosen moduli set. That is, the distance between each pair of elements must be derived as a 
function of the PRP modulus, m^, so that the corresponding amplitude response has die 
correct number of folds. For a given modulus mj, the number of folds needed to cover the 
dynamic range is given by 


•^i.SNS 


dynamic range M 

2 m; 2mi ’ 


(3-3) 


where M is given by Equation (3-2). Next consider the array factor given by Equation (2- 
21). Define ^=i)//2 given by 


^ = —sine. (3-4) 

From Equation (3-4) we have that: 


e=o' 

=> 

^=0 


o 

o 

11 

<x> 

=> 

7U ci 

(3-5) 

CD 

II 

o 


t d 



X ' 
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Each fold of the array factor corresponds to a change in phase of 2k which implies 

^(90°)-^(-90“) = 27c. (3-6) 

Substituting in Equation (3-6) the values of ^(±90°) found in Equation (3-5) we obtain 

2 jud ^ , . 

— = 2k, d = >.. (3-7) 

Therefore, the distance between elements is given by 

M 

2mj 

where X is the wavelength corresponding to the center frequency used to construct the 
antenna. For example, if an antenna is to be designed assuming a center frequency fQ=8.5 
GHz using mj=3, and 013 = 5 , the distance required between each pair of elements is: 

dj = 10 X, = 35.29 cm (mj = 3 ) 
dj = 7.5 X, = 26.47 cm {m^ = 4) 
dj = 6 X, =21.18 cm (m 3 = 5 ). 

With the correct distances between elements, the amplitude response from each pair of 
elements will be quantized using nij-l comparators. The function of the comparators is to 
generate the integers within each modulus by quantizing the amplitude response according 
to the modulus mj with which it is associated. Finally the recombination of the SNS code 
from each pair of elements will be the instantaneous direction of arrival. 
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C. SIMULATION RESULTS 


The array factor plot for the mj=3 array as the input direction of arrival varies from 
-90° to 90° is shown in Figure 3.3, As shown the folding period changes as the angle 6 
goes from 0° to ± 90°, 



Figure 3.3: Array Factor as a Function of 6 for a Two-Element Array Antenna 

With d=7.5 X. 


The SNS code cannot be applied to this waveform in a straight forward maimer. However, 
if the array factor is examined as a function of sin(e) as shown in Figure 3.4, the array 
factor corresponds to a symmetrically folded sinusoidal waveform (-1 < sin(0) < 1) with a 
constant period, thereby satisfying the principal condition for applying the SNS code. 

The chosen moduli have a dynamic range M=60, therefore the quantization width, 
considering the field of view (from 0° to ± 90°) that corresponds to the variation in sinO 


25 
























(from 0 to ± 1), will be equal to 2/60=0.0333. For a true symmetrically folded sinusoidal 
waveform as a function of 0, with constant period, the quantization width would be 
constant, independent of the direction of arrival, and equal to 3° (180760). This, however, 
is not the case. Comparing Figures 3.3 and 3.4 it is possible never the less to establish a 
relationship between the constant quantization width (when the array factor is plotted as a 
function of sin(6)) and the corresponding quantization width as a function of 0 (when the 
array factor is plotted as a function of 0). 



Figure 3.4: Array Factor as a Function of sin(0) for a Two-Element Array Antenna With 

d=7.5 X. 


As a function of 0 the quantization width depends on the direction of arrival of the 
incoming signal and is not uniform. The transition points from one quantization level to 
another quantization level occur at an angle 0j given by 
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0, =sin~^(x) 


(3-9) 


where x is given by 


x = -l,-l-h|-,-l + 4,...,l. (3-10) 

M M 

Table 3.2 shows the resulting transition points, 0^, and the respective quantization widths. 

The array factor for each pair of elements {(1,4), (1,3), and (1,2)} corresponding 
to modulus mj=3, m2=4, and m3=5 respectively are shown in Figures 3.5, 3.6, and 3.7. 
Also shown are the states within each modulus. The array factor is phased at -90° in aU 
cases so that the SNS code ‘000’ will correspond to our &st quantization level at -90° (see 
Table 3.1). To obtain this proper phase relationship between pairs, a phase shift equal to 7 t, 
JC/2, and tc was added to elements (1,2), (1,3), and (1, 4) respectively. 

The next step is to separately quantize the array factor waveform corresponding to 
each pair of elements using the SNS encoding procedure explained in Section A of this 
chapter. The number of threshold levels (or comparators) depends on the modulus nrij that 
is being considered and is given by nij-l as previously noted. The matching threshold 
values for the comparators are found by simulation. The transition points from one 
quantization level to another occur at the values of 0, shown in Table 3.2. 
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sin'^ X 

(degrees) 

Quantization 

Width 

sin'^ X 

e . 

(degrees) 

Quantization 

Width 

- 1.0000 

-90.0000 


0.0000 

0.0000 

1.9102 

- 0.9666 

-75.1649 

14.8351 

0.0333 

1.9102 

1.9123 

- 0.9333 

-68.9605 

6.2044 

0.0666 

3.8102 

1.9166 


-64.1581 

4.8025 

0.1000 

5.7392 

1.9231 

- 0.8666 

-60.0736 

4.0845 

0.1333 

7.6623 

1.9318 

- 0.8333 

-56.4427 

3.6309 

0.1666 

9.5942 

1.9429 

- 0.8000 

-53.1301 

3.3126 

0.2000 

11.5370 

1.9564 

- 0.7666 

-50.0555 

3.0746 

0.2333 

13.4934 

1.9726 

- 0.7333 

-47.1666 

2.8889 

0.2666 

15.4660 

1.9916 

- 0.7000 

-44.4270 

2.7396 

0.3000 

17.4576 

2.0136 

- 0.6666 

-41.8103 

2.6167 

0.3333 

19.4712 

2.0390 

- 0.6333 

-39.2965 

2.5138 

0.3666 

21.5102 

2.0680 

- 0.6000 

-36.8699 

2.4266 

0.4000 

23.5782 

2.1011 

- 0.5666 

-34.5181 

2.3518 

0.4333 

25.6793 

2.1389 

- 0.5333 

-32.2310 

2.2872 

0.4666 

27.8181 

2.1819 

- 0.5000 

-30.0000 

2.2310 

0.5000 

30.0000 

2.2310 

- 0.4666 

-27.8181 

2.1819 

0.5333 

32.2310 

2.2872 

- 0.4333 

-25.6793 

2.1389 

0.5666 

34.5181 

2.3518 

- 0.4000 

-23.5782 

2.1011 

0.6000 

36.8699 

2,4266 

- 0.3666 

-21.5102 

2.0680 

0.6333 

39.2965 

2.5138 

- 0.3333 

-19.4712 

2.0390 

0.6666 

41.8103 

2.6167 

- 0.3000 

-17.4576 

2.0136 

0.7000 

44,4270 

2.7396 

- 0.2666 

-15.4660 

1.9916 

0.7333 

47.1666 

2.8889 

- 0.2333 

-13.4934 

1.9726 

0.7666 

50.0555 

3.0746 

- 0.2000 

-11.5370 

1.9564 

0.8000 

53.1301 

3.3126 

- 0.1666 

-9.5942 

1.9429 

0.8333 

56.4427 

3.6309 

- 0.1333 

-7.6623 

1.9318 

0.8666 

60.0736 

4.0845 

- 0.1000 

-5.7392 

1.9231 

0.9000 

64.1581 

4.8025 

- 0.0666 

-3.8102 

1.9166 

0.9333 

68.9605 

6.2044 

- 0.0333 

-1.9102 

1.9123 

0.9666 

75.1640 

14.8351 

mm 

0.0000 

1.9102 

1.0000 

90.0000 



Table 3.2: Equivalent Quantization Levels Transition Points and Quantization 

Level Width. 
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Figure 3.5: Array Factor Waveform and Respective SNS Quantization Code for the Pair of 
Elements (1,4) Corresponding to Modulus mj=3. 




Figure 3.6: Array Factor Waveform and Respective SNS Quantization Code for the Pair of 
Elements (1,3) Corresponding to Modulus m 2 = 4 . 
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Figure 3.7: Array Factor Waveform and Respective SNS Quantization Code for the Pair of 
Elements (1, 2) Corresponding to Modulus m3=5. 


The normalized comparator matching threshold values for each pair of elements 
corresponding to modulus 3,4, and 5 are shown in Table 3.3. 


Elements 

Modulus 

T1 

T2 

T3 

T4 

(1, 2) 

5 

-0.8090 

-0.3090 

0.3090 

0.8090 

(1, 3) 

4 

-0.7071 

0.0000 

0.7071 


(1,4) 

3 

-0.5000 

0.5000 

1 

1 

1 

1 

1 

1 



Table 3.3: Normalized Comparator Matching Threshold Values for Each Pair of Elements. 
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The comparator states (or thermometer codes) from each channel are recombined in 
an SNS logic block to give the instantaneous resolved direction of arrival. Figure 3.8 
shovi's the resolved direction of arrival or antenna transfer function corresponding to the 60 
different possible combinations of the SNS for the moduli set (3, 4, 5). Note that the 
quantization width is not constant 



Figure 3.8: SNS Resolved Direction of Arrival for -90° < 6 < 90°(no element factor). 


The actual amplitude response of a linear array antenna is given by Equation (2-1) 
which means that a multiplicative factor (the element factor) has to be taken into 
consideration in addition to the array factor. The effect of the element factor is to introduce 
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a drop-off in the antenna gain. The element factor will depend on the type of elements that 
are to be used in the array antenna (waveguides, dipoles, etc.). For example consider that 
this element factor is characterized by a cosine waveform. Taking this characterization into 
account. Equation (2-1) can be rewritten as follows 

l^jcos(e). (3-11) 


E(6, ()) = 0) = Aq cos 


Figure 3.9 shows a plot of the electrical field response of the array assuming that the 
element factor is characterized by cos(0). The element factor modulates the array factor so 
that the resultant waveform no longer has a constant amplitude. 



Figure 3.9: Amplitude of the Electrical Field Response of a Two Linear Array: Elements 1 

and 2. 
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When the SNS quantization methodology is applied to the resultant waveform (same 
threshold levels) a number of errors result as shown in Figure 3.10. However, the results 
are still acceptable in a smaller interval (-20° < 6 < 20°). The antenna transfer function for 
this narrow region is shown in Figure 3.11. The encoding errors which exist at the code 
transition points are the result of some of the comparators not switching simultaneously as 
they should. These errors are easily isolated by using a few additional comparators in the 
minimum modulus channel [Ref. 5]. Investigation into the error isolation however, is 
beyond the scope of this thesis and is subject of future research. 



Figure 3.10: SNS Resolved Direction of Arrival for -90° < 0 < 90° (element factor 

included). 



























Figure 3.11: SNS Resolved Direction of Arrival for -20° < 0 < 20°(element factor 

included). 


D. EXPERIMENTAL RESULTS 

A linear antenna was built based on the theory developed previously using four 
rectangular waveguides mounted in a large ground plane. Figure 3.12 shows the 
configuration of the antenna (front plane, back plane, and mounted on the pedestal in the 
anechoic chamber). The connection between elements are as illustrated in Figure 3.2. 
Using the anechoic chamber facility at the Naval Postgraduate School, antenna pattern 
measurements (power response) for the H-field were taken for each pair of elements. 
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Figure 3.12.b: Photograph of the Built Antenna (Back Plane). 


35 














Figure 3.12.c: Photograph of the Built Antenna Mounted on the Pedestal in the 

Anechoic Chamber. 


Before evaluating the experimental results, the simulated results for this particular antenna 
(rectangular waveguide elements) are examined. The theoretical element factor (EF) for a 
rectangular waveguide when analyzing the H-field is given by 


gp _ cos(ka/2 sin6) 
Jt^ — (ka/2 sin0)^ 


COS0 


(3-12) 


where a is the width of the waveguide and k is the wave number [Ref. 3]. When Equation 
(3-12) is used as the element factor (instead of the cos(0) used previously in Section C), 
the resultant electrical field for each pair of elements is shown in Figure 3.13. Figure 3.14 
shows the resolved direction of arrival obtained when this electrical field is processed using 
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Electrical Field Response 


the SNS architecture presented in Section A of this chapter. The results compare favorably 
to those shown in Figure 3-10. 



Input Direction of Arrival 

Figure 3.13: Electrical Field Response for Mod 3, Mod 4, and Mod 5 When Using 

Rectangular Waveguide Elements. 
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Figure 3.14: SNS Resolved Direction of Arrival as a Function of 6, When Using the EF 

for a Rectangular Waveguide. 


The experimental measurements of the antenna pattern for each pair of elements is 
shown in Figure 3.15. Since the SNS architecture quantizes the amplitude of the received 
electrical field (not the power) it is necessary to convert to the collected measurements 
corresponding electrical field response. The relationship between the Power (P) and the 
Electrical Field (E) is given by 


P = 20 LogiolEl. 


(3-13) 
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Antenna Pattern (dB) 


The absolute value of the electrical field is then 


IEI = 10^'“. 


(3.14) 



Input Angle of Arrival 


Figure 3.15: Antenna Pattern for Each Pair of Elements (Experimental Data). 


Using Equation (3-14), the experimental measurements for each modulus (pair of elements) 
were used to calculate the absolute value of the electrical fields and are shown in Figure 
3.16. 
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Figure 3.17. Electrical field Response for Each Pair of Elements Corresponding to 

Modulus 3,4, and 5. 


The electrical field measurements for each pair of elements, representing modulus 3, 4, and 
5 is then quantized, and the output from each channel is recombined to obtain the resolved 
direction of arrival. The resolved direction of arrival is shown in Figure 3.18. 
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Figure 3.18: Resolved Direction of Arrival Using Experimental Data. 


Comparing the experimental results (Figure 3.18) with the simulated results (Figure 3.14) 
shows a good deal of agreement. However, considering just the interval from 0° to ± 20°, 
we notice many more encoding errors in the experimental results. If the electrical field was 
available directly in the anechoic chamber, the approximation would not have been 
necessary and the number of encoding errors would have been more in line with those of 
Figure 3.14. 


42 































IV. RESIDUE NUMBER SYSTEM ENCODING OF THE PHASE 

RESPONSE 

A. THE RESIDUE NUMBER SYSTEM (RNS) 

The ancient study of the Residue Number System begins with a verse from a third- 
century book, Suan-ching, by Sun Tzu [Ref. 7]: 

We have things of which we do not know the number. 

If we count them by threes, the remainder is 2, 

If we count them by fives, the remainder is 3, 

If we count them by sevens, the remainder is 2, 

How many things are there? 

The answer, 23. 

This 1700-year-old number system has been attracting a great deal of attention 
recently. Digital systems that use residue arithmetic units may play an important role in 
ultra-high speed, dedicated, real-time systems that support pure parallel processing of 
integer valued data. It is a “carry fi'ee” system that performs addition, subtraction, and 
multiplication as concurrent operations, side-stepping one of the principal arithmetic delays: 
managing carry information. Some applications such as fast Fourier transform, digital 
filtering, and image processing utilize the efficiencies of RNS arithmetic in addition and 
multiplication. The RNS has been employed efficiently in the implementation of digital 
processors. 

Sumlar to the SNS, the RNS is also composed of a number of pairwise relatively 
prime moduli rrij [Ref. 6]. The integers within each RNS modulus are representative of a 
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saw-tooth waveform with the period of the waveform equal to the modulus. In this 
scheme, given m, the integer values within the individual modulus are given by the row 
vector 


y = [0,l,...,m-l]. (4-1) 

Again, the required number of comparators for each channel is nij-l. Notice that all the 
integers within Equation (4-1) form a complete system of length 01 ^ (no ambiguities). For N 
equal to the number of moduli, the dynamic range of this scheme, as weU as the number of 
quantization levels, is given by 


M = nm,. (4-2) 

i=l 


Table 4.1 illustrates the RNS for the moduli mi=3, m2=4, m3=5. Figure 4.1 shows the 
RNS folding waveforms and the corresponding RNS states. 
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7 13 2 

8 2 0 3 

9 0 14 

10 1 2 0 

11 2 3 1 

12 0 0 2 


52 1 

53 2 

54 0 

55 1 

56 2 

57 0 

58 1 

5?_2 

60 0 



Table 4.1: An Example to Illustrate the RNS Encoding Procedure for 
the Moduli mi=3, m2=4 and m3=5. 
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□ : RNS states (threshold levels crossed) 



0 1 2345678 9 10 11 

Noimalized Input Direction of Arrival 

Figure 4.1: RNS Folding Waveforms and Output Codes for the Moduli 

m^=3, ni2=4, and m3=5. 

B, RNS ANTENNA ARCHITECTURE 

The phase response of a two-element interferometer corresponds to a saw-tooth 
folding waveform. The RNS antenna architecture being investigated is shown in Figure 
4.2. Phase detectors are used to combine the individual pairs of elements. Also, the 
comparator threshold levels and logic block are particular to the RNS (see Table 4.1). 
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DIRECTION OF ARRIVAL 0 


Figure 4.2: RNS Antenna Architecture. 


The distance between each pair of elements must be derived as a function of the 
PRP modulus, mj, so that the corresponding phase response has the correct number of 
folds. For a given modulus m^, the number of folds needed to apply the RNS code within 
the dynamic range is given by 
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^i.RNS ” 


_ dynamic range _ M 


m, 


m. 


(4-3) 


Next consider the two element interferometer shown in Figure 1.1. The phase difference 
between the received signal due to the separation between elements is given by 



A(j) = 

27cd . _ 

= —— sin6. 

X 

(4-4) 

From Equation (4-4) we have that: 




0 

= 0° 

=> A(j) = 0 


e 

= 90' 

., 27cd 

(4-5) 

0 

= -90° 

., 27td 

A^> = - . 



Since each fold of the phase response corresponds to a change in phase of 27t, we require 
that 


A(t)(90°) - A(l)(-90°) = 271 ( 4 - 6 ) 

which implies (from Equation (4-5) , 



Finally, the distance between elements is given by 
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(4-8) 



mi 



Using the same moduli set with center frequency of fo=8.5 GHz, the distances for each pair 
of elements are 


di =10X, 

= 7.5 X,, 

dj = 6.0 X. 

which are the same as those given for the SNS architecture. That is, the same hardware 
constructed for the SNS array can be used for the RNS array. 

C. SIMULATION RESULTS 

The correct distances between each pair of elements were found in Section B so that 
the output waveforms representative of modulus mi=3, m2=4, and m 3 = 5 , have the correct 
number of folds required to quantize the entire phase response within the field of view 
(from 0° to +90°). Figure 4.3 shows a plot of the phase response, as a function of the input 
direction of arrival for a two-element interferometer separated by a distance d=6X. The 
saw-tooth waveform stretches out as 0, the input direction of arrival, goes from 0° to ±90°. 
However, if the phase response is plotted as a function of sin(0), a perfect saw-tooth 
waveform is obtained as shown in Figure 4.4. 
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Phase Response (rad) Pti3se Response (rad) 



Figure 4.3; Phase Response as a Function of the Input Direction of Arrival for a 
Two-Element Interferometer with Separation d 3 = 6 X-. 



Figure 4.4: Phase Response as a Function of sin(0) for a Two-Element Interferometer with 

Separation d 3 = 6 X,. 
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Based on these results, the same analysis that was performed when dealing with the array 
factor can be applied here. In Figures 4.5, 4.6, and 4.7, the phase response waveform and 
respective output quantization code for modulus mj=3, modulus m2=4, and modulus m 3=5 
are plotted, as a function of the input direction of arrival. 




Figure 4.5: Phase Response Waveform and Respective RNS Quantization Code for the 
Pair of Elements (1,4) Corresponding to Modulus mj=3. 
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Figure 4.6; Phase Response Waveform and Respective RNS Quantization Code for the 
Pair of Elements (1,3) Corresponding to Modulus m 2 = 4 . 




Figure 4.7: Phase Response Waveform and Respective RNS Quantization Code for the 
Pair of Elements (1,2) Corresponding to Modulus m3=5. 
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The quantization of the phase response from each pair of elements is made 
separately using a small bank of comparators. The number of threshold levels (and 
comparators) are the same as that for the array factor amplitude analysis. However, the 
matching threshold values are different as shown in Table 4.2. 


Elements 

Modulus 

T1 

T2 

T3 

T4 

(1,2) 

5 

-1.8850 

-0.6283 

0.6283 

1.8850 

(1,3) 

4 

-1.5708 

0.0000 

1.5708 


(1,4) 

3 

-1.0472 

1.0472 




Table 4.2: Comparator Threshold Levels for Each Pair of Elements When 
Considering the Phase Analysis. 


For a specific input direction of arrival and pair of elements, there will be a 
corresponding integer (number of comparators in the on state) within the PRP modulus W;. 
The value of those integers depends on the measured phase difference associated with each 
pair of elements for that specific input direction of arrival. The comparator state for each 
pair together can then be represented by 3 RNS codes. The 3 RNS codes are indexed to an 
RNS logic block to give the resolved direction of arrival. Table 4.3 shows the 3 RNS 
codes and the resolved direction of arrival. 
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3RNS 

Codes 


000 


111 


222 


033 


104 


210 


021 


132 


203 


014 


120 


231 


002 


113 


224 


030 


101 


212 


023 


134 


200 


Oil 


122 


233 


004 


110 


221 


032 


103 


214 


Resolved DOA 
(degrees) 


10.5655 


31.1155 


33.3745 


40.5534 


45.7968 


62.1158 


3 RNS 

Resolved DOA 

Code 

(degrees) 

020 


131 

1 -72.0627 1 

202 


013 


124 

1 -58.2581 1 

230 

■K&Mi&lH 

001 

1 -51.5928 1 

112 

■BBsUiSH 

223 


034 


100 


211 

-38.0832 

022 

-35.6940 

133 


204 


010 


121 


232 

-24.6287 1 

003 


114 


220 

-18.4644 

031 

-16.4618 

102 

-14.4797 

213 


024 


130 


201 


012 


123 

-2.8664 1 

234 



Table 4.3: RNS Code and Equivalent Resolved Direction of Arrival for 
Modulus 3,4, and 5 Phased at 0° with Center Frequency fo=8.5 GHz. 
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Figure 4.8, shows a plot of the resolved direction of arrival as a function of 0, the 
input direction of arrival. This new system gives the direction of arrival information over 
the complete field of view without ambiguity. However the accuracy varies depending on 
the actual direction of arrival of the incoming signal. Notice also that the quantization width 
increases as the direction of arrival approaches ±90°. The quantization error wUl be 
addressed in Chapter V. 



—'-'-'-1-1-1_^_I_1_i_ 

-80 -60 -40 -20 0 20 40 60 80 

Input Direction of Arrival ( degrees ) 


Figure 4.8: Resolved Direction of Arrival in Function of the Input Direction of Arrival. 
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D. EXPERIMENTAL RESULTS 


Unfortunately it was not possible to obtain acceptable experimental results for the 
RNS antenna architecture. The main problem was the nonavailability of the phase 
detectors. Another problem was matching the elements with the correct length of cable in 
order to obtain the correct phase difference between each pair of elements. Future 
investigations are currently underway and will correct these problem. 
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V. FURTHER CONSIDERATIONS 


A. OVERVIEW 

Comparing the simulated results for the SNS antenna design with the results 
obtained for the RNS antenna design, it is clear that the second case permits acquiring the 
direction of arrival information over the complete field of view from 0° to ±90°. When 
using the SNS antenna it is impossible to obtain the direction of arrival information over the 
entire range from -90° to 90° due to the degradation in the amplitude response. Since the 
phase response has a “uniform” distribution (i.e., maximum and minimum with a constant 
value) over the complete field of view, it is possible to use the RNS code to quantize that 
signal correctly. Accordingly, only the phase analysis will be considered hereafter since the 
results shown are more encouraging. At this point we have not presented any results 
concerning the following aspects: 

• quantization errors associated with this new antenna design, 

• the antenna frequency response (incoming signal frequency different than the one 
used in the design), 


• the required size of the antenna related to the moduli set and the frequency used to 
build tile antenna, and 


• the influence of the chosen moduli set on the accuracy and quantization error of the 
antenna. 

These aspects are considered in this chapter for the RNS architecture (using the phase 
response). Some antennas design examples are given to show the degree of accuracy that 
can be reached with this direction finding architecture. 
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B. QUANTIZATION ERROR ASSOCIATED WITH THE RNS 
ARCHITECTURE 

The conversion of an analog (continuous) signal into a digital (discrete) 
representation is called quantization. Graphically, the quantization process means that a 
straight line representing the relation between the input and output of a linear continuous 
system is replaced by a staircase characteristic [Ref. 8], as shown in Figure 4.8. The 
difference of two adjacent discrete values is called the quantum or step size. Signals applied 
to a quantizer, with the input-output characteristic of Figure 4.8, are sorted into amplitude 
slices, and all input signals within plus or minus half a quantum step of the midvalue of a 
slice are replaced in the output by the midvalue. The quantizer error consists of the 
difference between the input and output signals of the quantizer so that the maximum 
instantaneous value of this error is half of one quantum step, and the total range of variation 
is from minus half a step to plus half a step. By analyzing Figure 4.8 again, it is clear that 
the quantization error for the RNS antenna is not constant since the code widths are not 
constant. The Appendix contains the Matlab™ program code to calculate the quantization 
error for a predetermined moduli set as a function of the input direction of arrival. Figure 
5.1 shows a plot of the characteristic quantization error curve associated with the resolved 
direction of arrival from Figure 4.8 for the moduli set (3,4,5) and for a frequency equal to 
8.5 GHz. For this case, the results show that the quantization error varies from about ±1° 
over the input range ±20°, to a maximum of ±7.42° for direction of arrivals near ±90°. 
Figure 5.2 it is a plot of both the resolved direction of arrival and the associated 
quantization error as a function of the input direction of arrival. 
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Figure 5.1: Characteristic Quantization Eiror Associated With the RNS Antenna When 
Considering the Moduli Set (3,4,5) and the f =fo=8.5 GHz. 



Figure 5.2: Resolved Direction of Arrival and Quantization Error. 
































C. RNS ANTENNA FREQUENCY RESPONSE 

A frequency (4) was chosen to calculate the required distance between each pair of 
elements so that the correct number of folds (for each pair of elements) were available to 
apply the RNS code and process (quantize) correctly the phase response without 
ambiguities over the entire field of view. The simulated results shown previously are 
obtained assuming that the frequency of the incoming signal (Q was equal to fp. In this 
case, the only error present is the quantization error. It is interesting to examine what 
happens when the frequency of the incoming signal has some other value different from fo. 
Figure 5.3, shows the resolved direction of arrival for an incoming signal with a frequency 
f =4.25 GHz. Comparing these results with the results in Figure 4.8, we notice the 
following: the total number of quantization levels is reduced by half and the resolved 
direction of arrival has large errors as the input direction of arrival goes from 0° to ±90°. 



Figure 5.3: Resolved Direction of Arrival for an Incoming Signal With a Frequency 

f,=4.25 GHz. 
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Alternatively, Figure 5.4 shows the resolved direction of arrival when the incoming signal 
frequency f,=17 GHz. In this case, we notice that the total number of quantization levels 
doubles when compared with the simulated results shown in Figure 4.8 (f,=fo). As a result, 
the resolved directions of arrival are ambiguous because the same 3 RNS codes correspond 
to more that one direction of arrival. In addition the overall error is very high. 



Figure 5.4: Resolved Direction of Arrival for an Incoming Signal With a Frequency 

f, = 17 GHz. 


The RNS antenna was designed assuming a moduli set (3,4,5) and a frequency (fp) equal 
to 8.5 GHz. For a given PRP modulus itij. Equation (4-3) gives the required number of 
folds to apply the RNS code within the dynamic range without ambiguities. It was found to 
be equal to 20 for modulus 3. In addition, the distance between each pair of elements is 
given by Equation (4-4). This distance is fixed when the antenna is built On other hand, in 
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the presence of a signal with a frequency other than f^, and from Equation (1-1), the phase 
difference between the received signal due to the separation between elements is given by 

A(l) = —sin 9. (5-1) 

Consider modulus 3, d=10X-o. Substituting that value in Equation (5-1) gives 

Ai|) = 25^sin6. (5-2) 

K 

To analyze Equation (5-2), three distinct cases need to be considered according to the value 
of the ratio 

• Casel: VX, = 1 =» f^ = f,. 

This can be considered the ideal case since the corresponding phase response (sawtooth 
waveform) will have the exact number of folds required to apply the RNS code within the 
dynamic range without ambiguities over the complete field of view (from 0° to ±90°). 

• Case 2; V^r<l =» fr<fo- 

In this case the phase response waveform will have less folds than the required number to 
apply the RNS code within the dynamic range. No ambiguities will be present; however, 
the quantization error will increase since a smaller number of quantization levels will exist 
to cover the same field of view (from 0° to ±90°). This case is directly related to Figure 5.3 
where the ratio X.oA.,=0.5. The equivalent number of folds and quantization levels are one 
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half those found for the ideal case. Figure 5.5 shows the phase response for both 
situations: f =8.5 GHz and f =4.25 GHz for modulus m3=5. 



Input Direction of Arrival ( degrees ) 


Figure 5.5: Phase Response Waveform (Mod 5) for f =8.5 GHz (Dashed) and 

f,=4.25 GHz. 


• Case 3: V\>1 

For this case the resultant phase response will have more folds than the required number to 
apply the RNS code. As a result we start to have ambiguous resolved directions of arrival 
as observed in Figure 5.4. Here, the ratio XJk=2', therefore the number of folds and 
quantization levels are twice the values found for the ideal case, as illustrated in Figure 5.6 
for modulus m3=5. 
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In all cases, the number of folds present in the phase response waveform will be 
proportional to the ratio 

In summary, for the results obtained in the third case, we won’t be able to find the 
direction of amval over the complete field of view without some ambiguous results. 
However, when the frequency of the incoming signal is smaller than fg, we do not have 
any ambiguities, but the overall error of the resolved direction of arrival increases. 



-80 -60 -40 -20 O 20 40 60 80 
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Figure 5.6: Phase Response Waveform (Mod 5) for f =8.5 GHz (Dashed) and 

f=17 GHz. 


This suggests that, for example, if one of the requisites of om antenna is to be able to find 
the direction of amval of any incoming signal with frequencies over some range, say 2 to 
18 GHz, the higher frequency in that range should be chosen to build our antenna, i.e., 18 
GHz for this case. By using this, the direction of arrival information will always be 
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obtained without ambiguities over the entire field of view. However, the overall 
quantization error will be different for each case depending on the number of quantization 
levels that will be used and, consequently, on the frequency value of the incoming signal. 

At this point, we determine how to correct the additional error when the frequency 
of the incoming signal is different than 4 (frequency used to build our antenna). To correct 
the error however, the firequency of the incoming signal must be known. With the 
measured frequency, the number of folds present in the phase response waveform is 
determined. The new number of folds present in the phase response over the entire field of 
view is given by 




ni.fn = 


\^rj 


n. 


i,RNS* 


(5-3) 


Consider modulus 3 and f^=4.25 GHz. For this particular case, the new number of folds 
(^,fn) to 10. From Equation (4-3) and knowing that the number of quantization 

levels is equal to the dynamic range, the actual number of quantization levels (Nqj) for a 
specific frequency is given by 


Nql= ni.f„ , (5-4) 

where n^^^ is the new number of folds and m^ represents the PRP modulus that is being 
considered. For our example, the number of quantization levels is equal to 30. With this 
number of quantization levels it is possible to find the transition points from one 
quantization level to another. The transition points will occur at some angle 0, for which the 
sin(0j) corresponds to some value in the range firom 0 to ±1 given by 
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(5-5) 


x(n) = ±- 

N 

■‘^QL 

for n = 0,1,2,..., Nql. Therefore, the transition angles are given by 

6,(n) = sin”^(x(n)). (5-6) 

Haying the transition angles, the midvalue of each quantization level can easily be found, as 
well as the corresponding quantization error. The new midvalues for each quantization level 
are then indexed to the RNS logic block so that for each 3 RNS codes will correspond the 
correct resolved direction of arrival. Since the phase response waveform is always phased 
at 0° (see Figures 4.5,4.6, and 4.7), it is possible to implement this new algorithm in our 
RNS logic block to obtain the correct resolved direction of arrival. Figure 5.7 shows the 
resolved direction of arrival and the quantization error for an incoming signal with a 
frequency f,=4.25 GHz, using this new algorithm to introduce the corrections. With this 
new algorithm, the error in the resolved direction of arrival when the frequency of the 
incoming signal is smaller than f^, becomes very small when compared with the results 
illustrated in Figure 5.3. Figure 5.8 shows the quantization error for this case using another 
scale. 


66 



Quantization Error (degrees) ^ Resolved DOA and Quantization Error (degrees 


Quantization Error 



-40 -20 o 20 40 

Input Direction of Arrival ( degrees ) 


igure 5.7: Resolved Direction of Arrival and Quantization Error for an Incoming Signal 

With a Frequency f^=4.25 GHz. 
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Figure 5.8: Characteristic Quantization Error Curve for an Incoming Signal With a 

Frequency f =4.25 GHz. 






























































D. RELATIONSHIP BETWEEN ANTENNA SIZE, MODULI 
SET AND FREQUENCY CHOSEN TO BUILD THE 
ANTENNA 

The size of the antenna is an important parameter when considering a new system to 
be implemented. The relationship between frequency and wavelength is given by X^c/f 
where c is the speed of light in free space. This relationship is important because the length 
of the DF antenna is given by the distance between the widest pair of elements which is 
proportional to the highest number of wavelengths (see Equation (4-8)) required to 
correctly apply the RNS code to a certain PRP moduli set. On other hand, when building 
the antenna, the highest frequency in the range of frequencies that are being observed must 
be chosen in order to be able to obtain the direction of arrival information over the entire 
field of view without ambiguities. Table 5.1 shows some different PRP moduli set 
combinations and the distance in wavelengths required for the widest pair of elements in the 
array. 


PRP Moduli Set 

(m„ mj 

Distance between the 
widest pair of elements 

(9, 13) 

6.5 X 

(17, 20) 

10 X 

(47, 50) 

25 X 

(97, 101) 

50.5 X 

(3, 4, 5) 

10 X 

(5, 7, 9) 

31.5 X 

(7, 9,11) 

49.5 X 

(11, 13, 15) 

97.5 X 


Table 5.1: Distance Between the Widest Pair of Elements 
According to the Moduli Set 
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Separation Between the Widest Pair of Elements (m) 


For each case considered in Table 5.1, Figure 5.9 shows the required distance between the 
widest pair of elements as a function of the frequency used to build the antenna. 



Figure 5.9: Required Distance Between the Widest Pair of Elements as a Function 
of the Frequency Used to Build the Antenna. 


The results show that it is advantageous to chose a frequency as high as possible to build 
the antenna. First, with a higher frequency it is possible to obtain the direction of arrival 
information without ambiguities over a larger bandwidth. In addition, for the same moduli 
set, the size of the antenna becomes smaller when using higher frequencies to build it. 
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E. QUANTIZATION ERROR BEHAVIOR VS MODULI SET 


Another interesting point is to investigate how the quantization error varies 
acx:ording to the chosen moduli set. The overall quantization error of the system is directly 
related to the chosen PRP moduli set, since it determines the dynamic range, and therefore 
the number of quantization levels available to quantize the phase response over the entire 
field of view. With a higher number of quantization levels, the quantization error is smaller. 
Therefore, a better performance of our antenna is to be expected. Figiures 5.10 through 
5.17 show the characteristic quantization error curve associated with each of the moduli set 
considered in Table 5.1, assuming that the antenna is build using a frequency (fo) equal to 
18 GHz. 



Input Direction of Arrival ( degrees ) 

Figure 5.10: Characteristic Quantization Error Curve for the PRP Moduli Set (9,13). 
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Figure 5.11: Characteristic Quantization Error Curve for the PRP Moduli Set (17,20). 
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Figure 5.12; Characteristic Quantization Error Curve for the PRP Moduli Set (47,50). 
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Figure 5.13: Characteristic Quantization Error Curve for the PRP Moduli Set (97,101). 
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Figure 5.15: Characteristic Quantization Error Curve for the PRP Moduli Set (5,7, 9). 
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Figure 5.17: Characteristic Quantization Error Curve for the PRP Moduli Set (11,13,15). 


Comparing the results obtained for each case it is clear that when using a N moduli 
set with N=2, the quanti25ation error is smaller than when using a moduli set with N=3, 
assuming that the widest distance between elements is the same. Therefore, with a moduli 
set with N=2, it is possible to obtain a very high accuracy with relative small antenna sizes. 
However, the required number of comparators for this scheme (N=2) becomes very high 
when compared with the scheme using a moduli set with N=3. Table 5.2 contains a 
compilation of aU information of each case considered in Table 5.1. In addition, the last 
column shows the equivalent rms accuracy of an interferometer (conventional methods) 
assuming that the direction of arrival is equal to 0° and the SNR is equal to 0 dB. The rms 
angular accuracy of an interferometer in radians [Ref. 1] is given by 


_ Ag 

TtVsr^’ 


(5-7) 
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where Aa is the separation between adjacent nulls given by 


Aa = 


X 

dcosG 


(5-8) 


Moduli Set 

Largest 

dist. 

— w. — 

Comparators 

Min. Quant. 
Error (deg.) 
f,= f,=18 GHz 

Min. Quant. 
Error (deg.) 
f, = 2 GHz 


(9,13) 

6.5 X 

20 

± 0.4897 

± 4.425 

4.897 

(17,20) 

10 X 

35 

±0.1685 

±1.517 

1.827 

(47,50) 

25 X 

95 

± 0.0244 

± 0.2194 

0.729 

(97,101) 

50.5 X 

196 

± 0.0058 

± 0.0526 

0.365 

(3,4,5) 

10 X 

9 

± 0.9551 

± 8.729 

1.827 

Ksin 

31.5 X 

18 

±0.1819 

±1.638 

0.579 

(7,9,11) 

49.5 X 

24 

± 0.0827 

± 0.744 

0.368 

(11,13,15) 

91.5 X 

36 

± 0.0267 

± 0.2404 

0.187 


Table 5.2; Comparison Between Different RNS Architecture Schemes and the 

Conventional Methods. 


As shown in Table 5.2, with this new architecture it is possible to obtain the direction of 
amval information with a substantial improvement in accuracy when compared with the 
conventional interferometer methods. 
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VI. SUMMARY AND CONCLUDING REMARKS 


In this thesis two novel architectures for high-resolution wideband direction finding 
based in the SNS and RNS were developed and investigated. The first architecture 
incorporates a SNS encoding of the amplitude of the electrical field response (symmetrical 
folding waveform) of the phase response (sawtooth folding waveform) array. The second 
architecture incorporates a RNS encoding of the interferometer phase. In addition, an array 
antenna was constmcted in a ground plane using rectangular waveguide elements with a 
center frequency f(,=8.5 GHz and the SNS architecture was tested. 

The simulated results show that the SNS architecture is sensitive to the drop off in 
the antenna gain due to the element factor effect associated with the kind of individual 
elements used to build the antenna. In this case it is not possible to obtain the direction of 
arrival information over the entire field of view (fi'om 0° to ±90°). However, we still have a 
small interval (from 0° to ±20°) where it is possible to have that information. The 
experimental results indicate a considerable agreement for this case. 

For the RNS architecture the simulated results indicate that it provides a more 
robust solution to high-resolution DF having an excellent performance. In fact, since the 
phase response has an uniform” distribution it is possible to obtain the direction of arrival 
information over the entire field of view. To avoid ambiguities in the resolved direction of 
arrival (over the entire field of view) the spacing between elements must be established for 
the highest frequency to be received. On other hand, the frequency of the incoming signal 
has to be known in order to make the necessary corrections in the RNS logic block using a 
fast correction algorithm to obtain the correct direction of arrival. The performance of the 
algorithm was quantified and shown to have excellent results over a large bandwidth. 
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The quantization error associated with both architectures depends on the moduli set 
that is chosen to build the antenna, which also determines the size of the antenna and the 
required number of comparators. Table 5.2, in Chapter V, illustrates some specifications 
for the RNS architecture using different moduli set combinations. It also shows the 
equivalent rms accuracy for the conventional methods according to the distance between the 
widest pair of elements. Comparing those values with the values found for the RNS 
architecture we notice a great improvement in accuracy with the new architecture. 

The advantages of the RNS architecture over the conventional methods can be 
summarized as follows: 

• Higher resolution, higher bandwidth capability, 

• Small number of antenna elements, 

• Small array size to achieve higher accuracy, 

• Only few comparators are required, depending on the chosen moduli set, 

• Full coverage firom 0° to ±90° without ambiguities. 

It was not possible to obtain acceptable experimental results for the RNS antenna 
architecture. The main problem was the nonavailability of the phase detectors. However, 
future investigations are currently underway and will correct this and another problems. 
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APPENDIX: MATLAB’^^ CODE 


%Program 1 

%This program 'Ampanal.m' calculates the resolved Direction of Arrival using the SNS 
%code to process the amplitude response of the array for the moduli set (3,4,5). 

clear aU 
whitebg 
global ml 
global m2 
global m3 
global fo 
global ff 
global ds 
global rad 
global M 
global theta 

obj=input(’Do you want the Array Factor (1) or the Electrical field (2) analysis ?’) 
ds=input('What is the step size? ') 

fo=inputCWhat was the frequency used to built the antenna (GHz) ?') 
fr=input('What is the frequency of the received signal (GHz) ?’) 


ml =3; 

m2=4; 

m3=5; 

fo=fo*le9; 

fr=fr*le9; 


rad=pi/180; 

theta=-90:ds:90; 

c=3e8; 

L=c/fo; 

Ll=c/fr; 


M=ml*m2*m3; 


^**********=ic********£)jgj^j^j>gg Between Each Pair of Elernents****************** 

dl=0.5*M/ml*L; 

d2=0.5*M/m2*L; 

d3=0.5*M/m3*L; 


%*****Generates the Array Factor or the Electrical field for Each Pair of Elements****** 

phil=(2*pi*dl/L)*sin(theta*rad); 

phi2=(2*pi*d2/L)*sin(theta*rad); 

phi3=(2*pi*d3/L)* sin(theta*rad); 


if obj==l 

modlA=cos(0.5*(phil)+pi); 
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mod2A=cos(0.5*(phi2)+pi/2); 

mod3A=cos(0.5*(phi3)+pi); 

else 

modlA=cos(0.5*(phil)+pi) *cos(theta*rad); 
mod2A=cos(0.5*(phi2)+pi/2).*cos(theta*rad); 
mod3A=cos(0.5*(phi3)+pi).*cos(theta*rad); 
end 

figure(l) 
subplot(311) 
plot(theta,mod 1 A) 
axis([-90,90,-l,l]); 

xlabel('Input Direction of Arrival (degrees)') 
grid 

subplot(312) 

plot(theta,mod2A) 

axis([-90,90,-l,l]); 

xlabel('Input Direction of Arrival (degrees)') 

ylabelCArray Factor Amplitude') 

grid 

subplot(313) 

plot(theta,mod3A) 

axis([-90,90,-l,l]); 

xlabel(’Input Direction of Arrival (degrees)') 
grid 

^******;i!*****!i!*(^^Qyj^^g Threshold Values for Each Modulus**************** 
thresa 

^*********( 3 gjjgj.a(gg {he Comparators RNS Output code for Modulus nil=3********** 
for i=l:180/ds+l 
ifmodlA(i)<=:Tl(l) 
modl(i)=0; 

elseif modlA(i)>Tl(2) 

modl(i)=2; 

else 

modl(i)=l; 

end 

end 

end 

%*********Generates the Comparators RNS Outout code for Modulus nT2=4********** 

for i=l:180/ds+l 

ifmod2A(i)<=T2(l) 

mod2(i)=0; 

elseif mod2A(i)<=T2(2) 
mod2(i)=l; 

elseif mod2A(i)<=T2(3) 
mod2(i)=2; 
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else 

mod2(i)=3; 

end 

end 

%*********Generates the Comparators RNS Output code for Modulus m 3 = 5 ********** 

for i=l:180/ds+l 

ifmod3A(i)<=T3(l) 

mod3(i)=0; 

elseifmod3A(i)<=T3(2) 

mod3(i)=l; 

elseifmod3A(i)<=T3(3) 

mod3(i)=2; 

elseif mod3A(i)<=T3(4) 

mod3(i)=3; 

else 

mod3(i)=4; 

end 

end 

%******Calculates the Resolved Direction of Arrival for the moduli Set (3,4, 5 )******** 
SNSLogicBlock 

%*Plots the Resolved Direction of Arrival as a Function of the Input Direction of Arrival** 

figure( 2 ) 

plot(theta,th,’g') 

xlabel('Input Direction of Arrival (degrees)') 
ylabel('Resolved Direction of Arrival (degrees)') 
grid 

axis([-90 90 -90 90]) 
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% Program 2 

% This program ‘thresa.m’ calculates the Threshold Values for the moduli set (3, 4, 5) 
%when analyzing the Amplitude Response. 

global ml 
global m2 
global m3 
global fo 
global ds 
global M 
global L 


%*******Calculates the transition angles from one to another quantization level ******** 

x=-l:2/M:l; 

y=asin(x); 


^*************5t:Qj^jQjjj^^gg Threshold Values for modulus nil=3*************** 

for l=l:ml-r 

angl(l)=y(l+l); 

phi 1 (l)=(2*pi*d 1/L) *sin(ang 1 (1)); 

Tl(l)=cos(0.5*(phil(l))+pi); 

end 

^^**************(^^g^^^g 5 Threshold Values for Modulus ni2=4*************** 

forr=l:m2-l 

ang2(r)=y(r+l); 

phi2(r)=(2*pi*d2/L)*sin(ang2(r)); 

T2(r)=cos(0.5*(phi2(r))+pi/2); 

end 

%**************calculates the Threshold Values for Modulus m3=5*************** 

for h=l:m3-l 

ang3(h)=y(h+l); 

phi3(h)=(2*pi*d3/L)’^sin(ang3(h)); 

T3(h)=cos(0.5*(phi3(h))+pi); 

end 
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%Program 3 

%This program ‘SNSLogicBlock.m’ calculates the resolved direction of arrival according 
to the 3-SNS code. 


global ds 
global ml 
global m2 
global m3 
global fo 
global fr 
global M 
global theta 

%****Combines the SNS Code coming from each channel ml=3, m2=4, and m3=5***** 
res=[modr mod2' mod3’]; 

%***Calculates the Resolved Direction of Arrival according to the 3-SNS Code value**** 
for u=l:180/ds+l 
if res(u,:)==([0,0,0]) 
th(u)=-82.5824; 
elseif res(u,:)==([l,l,l]) 
th(u)=-72.0627; 
elseif res(u,:)=([2,2,2]) 
th(u)=-66.5593; 
elseif res(u,:)=([2,3,3]) 
th(u)=-62.1158; 
elseif res(u,:)=([l,3,4]) 
th(u)=-58.2581; 
elseif res(u,:)==([0,2,4]) 
th(u)=-54.7864; 
elseif res(u,:)==([0,l,3]) 
th(u)=-51.5928; 
elseif res(u,:)==([!,0,2]) 
th(u)=-48.6110; 
elseif res(u,;)==([2,0,l]) 
th(u)=-45.7968; 
elseif res(u,:)==([2,l,0]) 
th(u)=-43.1187; 
elseif res(u,:)=([!,2,0]) 
th(u)=-40.5534; 
elseif res(u,:)=([0,3,l]) 
th(u)=-38.0832; 
elseif res(u,:)==([0,3,2]) 
th(u)=-35.6940; 
elseif res(u,:)=([l,2,3]) 
th(u)=-33.3745; 
elseif res(u,:)~([2,l,4]) 
th(u)=-31.1155; 
elseif res(u,:)=([2,0,4]) 
th(u)=-28.9091; 
elseif res(u,:)==([l,0,3]) 
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th(u)=-26.7487; 

elseif res(u,:)=([0,l,2]) 

th(u)=-24.6287; 

elseif res(u,:)=([0,2,l]) 

th(u)=-22.5442; 

elseif res(u,:)==([!,3,0]) 

th(u)=-20.4907; 

elseif res(u,:)==([2,3,0]) 

th(u)=-18.4644; 

elseif res(u,:)==([2,2,l]) 

th(u)=-16.4618; 

elseif res(u,:)=([l,l,2]) 

th(u)=-14.4797; 

elseif res(u,:)=([0,0,3]) 

th(u)=-12.5152; 

elseif res(u,:)==([0,0,4]) 

th(u)=-10.5655; 

elseif res(u,:)=([!, 1,4]) 

th(u)=-8.6382; 

elseif res(u,:)=([2,2,3]) 

th(u)=-6.7007; 

elseif res(u,:)==([2,3,2]) 

th(u)=-4.7809; 

elseif res(u,;)=([l,3,l]) 

th(u)=-2.8664; 

elseif res(u,;)==([0,2,0]) 

th(u)=-0.9551; 

elseif res(u,:)==([0,l,0]) 

th(u)=0.9551; 

elseif res(u,:)==([!,0,1]) 

th(u)=2.8664; 

elseif res(u,:)=([2,0,2]) 

th(u)=4.7809; 

elseif res(u,;)==([2,l,3]) 

th(u)=6.7007; 

elseif res(u,:)==([ 1,2,4]) 

th(u)=8.6382; 

elseif res(u,;)==([0,3,4]) 

th(u)=10.5655; 

elseif res(u,:)==([0,3,3]) 

th(u)=12.5152; 

elseif res(u,;)==([l,2,2]) 

th(u)=14.4797; 

elseif res(u,:)==([2,l,l]) 

th(u)=16.4618; 

elseif res(u,:)==([2,0,0]) 


th(u)=l 8.4644; 


elseif res(u,;)=([ 1,0,0]) 
th(u)=20.4907; 


elseif res(u,:)==([0,l,l]) 
1ii(u)=22.5442; 


! 
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elseif res(u,:)==([0,2,2]) 

th(u)=24.6287; 

elseif res(u,:)=([!,3,3]) 

th(u)=26.7487; 

elseif res(u,:)=([2,3,4]) 

th(u)=28.9091; 

elseif res(u,:)==([2,2,4]) 

th(u)=31.1155; 

elseif res(u,:)==([!,1,3]) 

th(u)=33.3745; 

elseif res(u,:)=([0,0,2]) 

th(u)=35.6940; 

elseif res(u,:)==([0,0,l]) 

th(u)=38.0832; 

elseif res(u,:)==([l,l,0]) 

th(u)=40.5534; 

elseif res(u,;)==([2,2,0]) 

th(u)=43.1187; 

elseif res(u,:)==([2,3,l]) 

th(u)=45.7968; 

elseif res(u,:)==([l,3,2]) 

th(u)=48.6110; 

elseif res(u,:)==([0,2,3]) 

th(u)=51.5928; 

elseif res(u,:)=([0,l,4]) 

th(u)=54.8;%54; 

elseif res(u,:)=([l,0,4]) 

th(u)=58.2581; 

elseif res(u,:)==([2,0,3]) 

th(u)=62.1158; 

elseif res(u,:)==([2,l,2]) 

th(u)=66.5593; 

elseif res(u,:)=([l,2,l]) 

th(u)=72.0627; 

elseif res(u,:)=([0,3,0]) 

th(u)=82.5824; 

end 

end 
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%Program 4 

%This program 'Phanalnew' calculates the Resolved Direction of Arrival using the RNS 
%code to process the phase response of the array for the moduli set (3,4, 5). 

clear all 
whitebg 
global ml 
global m2 
global m3 
global fo 
global fr 
global ds 
global rad 
global M 
global theta 

ds=input('What is the step size? ') 

fo=input('What was the frequency used to built the antenna (GHz) ?') 
fr=input(’What is the frequency of the received signal (GHz) ?') 


ml =3; 

m2=4; 

m3=5; 

fo=fo*le9; 

fr=fr*le9; 

rad=pi/180; 

theta=-90;ds:90; 

c=3e8; 

L=c/fo; 

Ll=c/fr; 


M=ml*m2*m3; 


%*************^=*****£)jj,^gjj^,gg Between Each Pair of Elernents****************** 

dl=0.5*M/ml*L; 

d2=0.5*M/m2*L; 

d3=0.5*M/m3*L; 


^**************Generates the Array Factor for Each Pair of Elements************** 

phil=(2*pi*dl/Ll)*sm(theta*rad); 

modlp=angle(exp(j*((phil+pi)))); 

phil=(2*pi*d2/Ll)*sin(theta*rad); 
mod2p=angle(exp(j *((phi 1 )+pi))); 

phi I=(2*pi*d3/Ll )*sin(theta*rad); 
mod3p=angle(exp(j *((phi 1 )+pi))); 
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figure(l) 

subplot(311) 

plot(theta,modlp) 

axis([-90,90,-4,4]); 

xlabel('Input Direction of Arrival (degrees)') 
grid 

subplot(312) 

plot(theta,mod 2 p) 

axis([-90,90,-4,4]); 

xlabel('Input Direction of Arrival (degrees)') 

ylabel('Phase Response (pi)') 

grid 

subplot(313) 

plot(theta,mod3p) 

axis([-90,90,-4,4]); 

xlabel('Input Direction of Arrival (degrees)') 
grid 

<^************** 0 ^jj.yj^^gg Threshold Values for Each Modulus**************** 
thres 

%*********Generates the Comparators RNS Output code for Modulus ml= 3 ********** 
for i=l:180/ds+l 
ifmodlp(i)<=Tl(l) 
modl(i)= 0 ; 
elseif modlp(i)>Tl(2) 
modl(i)= 2 ; 
else 

modl(i)=l; 

end 

end 

end 

%*********Generates the Comparators RNS Output code for Modulus m 2 = 4 ********** 

for i=l:180/ds+l 

ifmod2p(i)<=T2(l) 

mod 2 (i)= 0 ; 

elseif mod2p(i)<=T2(2) 
mod 2 (i)=l; 

elseif mod2p(i)<=T2(3) 

mod 2 (i)= 2 ; 

else 

mod2(i)=3; 

end 

end 

%*********Generates the Comparators RNS Output code for Modulus m 3 = 5 ********** 
for i=l:180/ds+l 
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if mod3p(i)<=T3(l) 
mod3(i)=0; 

elseif mod3p(i)<=T3(2) 
mod3(i)=l; 

elseif mod3p(i)<=T3(3) 
mod3(i)=2; 

elseif mod3p(i)<=T3(4) 

mod3(i)=3; 

else 

mod3(i)=4; 

end 

end 

%******Calculate the Resolved Direction of Arrival and the Quantization Error********* 
^************************fQj. rnoduli Set (3,4, S)**********************!):** 
RNSlogicBlock 

%*********Plots the Resolved Direction of Arrival and the Quantization Error********* 

^*****************^5 ^ function of the Input Direction of Arrival***************** 

figure(2) 

plot(theta,th,'g') 

hold 

xlabel(’Input Direction of Arrival (degrees)') 

ylabel('Resolved DOA and Quantization Error (degrees)') 

axis([-90 90 -90 90]) 

plot(y([2:size(y,2)]),k,’r’) 

grid 
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%Program 5 

%This program ‘thres.m’ calculates the Threshold values for the moduli set (3, 4, 5) 
%when analyzing the phase response. 

global ml 
global m2 
global m3 
global fo 
global fr 
global ds 
global M 
global L 


%*******Calculate the transition angles from one to another quantization level********* 

x=-l:2/M:l; 

y=asin(x); 


%:ic*************cal(;.uiate the Threshold Values for modulus nil=3**************** 

for l=l:ml-l 

angl(l)=y(l+l); 

phi 1 (l)=(2*pi*dl/L)*sin(ang 1 (1)); 

T1 (l)=angle(exp(j*((phi 1 (l))+pi))); 
end 

^**************Caiculate the Threshold Values for modulus ni2=4**************** 

for r=l:m2-l 

ang2(r)=y(r+l); 

phi2(r)=(2*pi*d2/L)*sin(ang2(r)); 

T2(r)=angle(exp(j*((phi2(r))+pi))); 

end 


^******<!*******QQculate the Threshold Values for modulus ni3=5**************** 

forh=l:m3-l 

ang3(h)=y(h+l); 

phi3(h)=(2*pi*d3/L)*sin(ang3(h)); 

T3(h)=angle(exp(j*((phi3(h))+pi))); 

end 
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% Program 6 

%This pro^am ‘RNSlogicBlock.m’ calculates the Resolved Direction of Arrival and 
%Quantization Error according to the 3-RNS code. It also uses an algorithm to calculate the 
%correct quantization level midvalue when the frequency of the received signal (fr) is 
%different form the center frequency (fo) used in the antenna design. 


global rad 
global ml 
global m2 
global m3 
global fo 
global ff 
global M 
global theta 


%*****Calculate the total number of quantization levels over the entire field of view***** 

%*******=i:***and the transition angles from one to another quantization level********** 

nf=fr/fo*M/ml; 

nql=nf*ml; 

p=0:2/nql:l; 

if max(p)~=l 

p=[p 11; 

end 

p=[ -rot90(p(2:size(p,2)))’ p]; 
trang= l/rad*asin(p); 

%**Calculate the quantization level width and the midvalues for each quantization level*** 
for n=l:size(p,2)-l 

midval(n)=trang(n)+(trang(n+l)-trang(n))/2; 

qlw(n)=midval(n)-trang(n); 

end 

if rem(size(midval,2),2)==0 

a=midval; 

b=size(a,2)/2; 

end 

if rem(size(midval,2),2)==l 

a=midval; 

b=(size(a,2)+l)/2; 

end 

if size(a,2)~=60 
a=[zeros(l,b) a zeros(l,b)]; 
end 

b=size(a,2)/2; 

for n=2:2*size(trang,2)-l 

if rem(n,2)==0 

y(n)=trang((n)/2); 

else 

y(n)=trang((n+l)/2); 

end 
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end 


^*********x:***********%Qg2cuiate the Quantization Error*********************** 

for n=l:size(trang,2)-l 

c(n)=(trang(n+l)-trang(n))/2; 

end 

for r=l:2*size(trang,2)-2 
if rem(r,2)==0 
k(r)=c(r/2); 

k(r)=-c((r+l)/2); 

end 

end 


%****Combmes the RNS Code corning from each channel ml=3, m2=4, and m3=5***** 

res=[mod3' mod4' mod5’]; 

for u=l:180/ds+l 
ifres(u,:)==([0,0,0]) 
th(u)=a(b+l); 
elseif res(u,:)==([!,1,1]) 
th(u)=a(b+2); 
elseif res(u,:)=([2,2,2]) 
th(u)=a(b+3); 
elseif res(u,:)~([0,3,3]) 
th(u)=a(b+4); 
elseif res(u,:)=([l,0,4]) 
th(u)=a(b+5); 
elseif res(u,:)==([2,1,0]) 
th(u)=a(b+6); 
elseif res(u,;)==([0,2,l]) 
th(u)=a(b+7); 
elseif res(u,:)==([l,3,2]) 
th(u)=a(b+8); 
elseif res(u,:)=([2,0,3]) 
th(u)=a(b+9); 
elseif res(u,:)=([0,l,4]) 
th(u)=a(b+10); 
elseif res(u,:)==([l,2,0]) 
th(u)=a(b+l 1); 
elseif res(u,:)==([2,3,l]) 
th(u)=a(b+12); 
elseif res(u,:)==([0,0,2]) 
th(u)=a(b+13); 
elseif res(u,:)==([l,l,3]) 
th(u)=a(b+14); 
elseif res(u,:)==([2,2,4]) 
th(u)=a(b+15); 
elseif res(u,:)=([0,3,0]) 
th(u)=a(b+16); 
elseif res(u,:)=([l,0,l]) 
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th(u)=a(b+17); 

elseif res(u,:)==([2,l,2]) 

th(u)=a(b+18); 

elseif res(u,;)==([0,2,3]) 

th(u)=a(b+19); 

elseif res(u,:)==([!,3,4]) 

tii(u)=a(b+20); 

elseif res(u,:)==([2,0,0]) 

th(u)=a(b+21); 

elseif res(u,:)==([0,l,l]) 

th(u)=a(b+22); 

elseif res(u,:)==([!,2,2]) 

th(u)=a(b+23); 

elseif res(u,:)==([2,3,3]) 

th(u)=a(b+24); 

elseif res(u,:)==([0,0,4]) 

th(u)=a(b+25); 

elseif res(u,:)==([!, 1,0]) 

th(u)=a(b+26); 

elseif res(u,:)==([2,2,l]) 

th(u)=a(b+27); 

elseif res(u,:)==([0,3,2]) 

th(u)=a(b+28); 

elseif res(u,:)==([l,0,3]) 

th(u)=a(b+29); 

elseif res(u,:)==([2,l,4]) 

th(u)=a(b+30); 

elseif res(u,:)==([0,2,0]) 

th(u)=a(b-29); 

elseif res(u,:)==([l,3,l]) 

th(u)=a(b-28); 

elseif res(u,:)==([2,0,2]) 

th(u)=a(b-27); 

elseif res(u,:)==([0,l,3]) 

th(u)=a(b-26); 

elseif res(u,:)==([l,2,4]) 


th(u)=a(b-25); 

elseif res(u,:)==([2,3,0]) 

th(u)=a(b-24); 


elseif res(u,:)==([0,0,l]) 
th(u)=a(b-23); 


elseif res(u,:)==([!,1,2]) 
th(u)=a(b-22); 


elseif res(u,:)==([2,2,3]) 
th(u)=a(b-21); 


elseif res(u,:)==([0,3,4]) 
th(u)=a(b-20); 
elseif res(u,:)==([l,0,0]) 
th(u)=a(b-19); 


elseif res(u,:)==([2,l,l]) 
th(u)=a(b-18); 



elseif res(u,:)=([0,2,2]) 

th(u)=a(b-17); 

elseif res(u,:)==([l,3,3]) 

th(u)=a(b-16); 

elseif res(u,:)==([2,0,4]) 

th(u)=a(b-15); 

elseif res(u,:)==([0,l,0]) 

th(u)=a(b-14); 

elseif res(u,:)==([l,2,l]) 

th(u)=a(b-13); 

elseif res(u,:)=([2,3,2]) 

th(u)=a(b-12); 

elseif res(u,:)==([0,0,3]) 

th(u)=a(b-ll); 
elseif res(u,:)==([!,1,4]) 
th(u)=a(b-10); 
elseif res(u,:)==([2,2,0]) 
th(u)=a(b-9); 
elseif res(u,:)==([0,3,l]) 
th(u)=a(b-8); 
elseif res(u,;)==([l,0,2]) 
th(u)=a(b-7); 
elseif res(u,:)==([2,l,3]) 
th(u)=a(b-6); 
elseif res(u,;)==([0,2,4]) 
th(u)=a(b-5); 
elseif res(u,:)==([!,3,0]) 
th(u)=a(b-4); 
elseif res(u,:)==([2,0,l]) 
th(u)=a(b-3); 
elseif res(u,:)=([0,l,2]) 
th(u)=a(b-2); 
elseif res(u,:)=([l,2,3]) 
th(u)=a(b-iy, 
elseif res(u,:)=([2,3,4]) 
th(u)=a(b); 
end 
end 
end 
end 




%Program 7 

%This program 'Quanterror' calculates the Quantization Error for any PRP N moduli set, 
%where N is equal to 2 or 3, and any frequency. 

clear aU 
whitebg 

ds=input('What is the step size? ') 

N=input('What is the moduli set size (2 or 3) ?') 
if N==2 

ml=input('What is the value of modulus ml ?') 
m2=input('What is the value of modulus m2 ?') 
end 

ifN=:=3 

ml=input('What is the value of modulus ml ?') 
m2=input('What is the value of modulus m2 ?’) 
m3=input('What is the value of modulus m3 ?') 
end 

fo=input('What was the frequency used to built the antenna (GHz) ?') 
fr=input(’What is the frequency of the received signal (GHz) ?') 

%**********Number of Folds Over the Entire field of View for Modulus ml********** 
nf=fr/fo*M/ml; 

%********Total Number of Quantization Levels Over the Entire Field of View********* 
nql=nf*ml; 


%************x=Gejierates the Quantization Level Transistion Angles**************** 

p=0:2/nql:l; 

if max(p)~=l 

p=[p 1]; 

end 

p=[ -rot90(p(2:size(p,2)))’ p]; 
trang= l/rad*asin(p); 

^***********si:**QjQj.jjjj^^gg Midvalue of Each Quantization Level**************** 
for n=l:size(p,2)-l 

midval(n)=trang(n)+(trang(n+l)-trang(n))/2; 

qlw(n)=midval(n)-trang(n); 

end 

^***********si!********(]; 3 jgy 2 jj^gg jjjg Quantization Error************************ 

for n=2:2*size(trang,2)-l 

if rem(n,2)==0 

y(n)=trang((n)/2); 

else 

y(n)=trang((n+l)/2); 

end 

end 
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for n=l:size(trang,2)-l 

c(n)=(trang(n+l)-trang(n))/2; 

end 

for r=l:2*size(trang,2)-2 
if rem(r,2)==0 
k(r)=c(r/2); 

k(r)=-c((r+l)/2); 

end 

end 

%**********piots thg Resolved Direction of arrival and Quantization Error************ 
^****************^5 Function of the Input Direction of Arrival****************** 
figure(l) 

plot(y([2:size(y,2)]),k,'r') 

xlabel('Input Direction of Arrival (degrees)') 

ylabelC Quantization Error (degrees)') 

axis([-90 90 -round(max(k)+0.5) round(max(k)+0.5)]) 

grid 
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